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Abstract 

Decentralized  control  methods  are  appealing  in  coordination  of  multiple  vehicles  due  to  their  low 
demand  for  long-range  communication  and  their  robustness  to  single-point  failures.  An  important 
approach  in  decentralized  multi- vehicle  control  involves  artificial  potentials  or  digital  pheromones.  In 
this  paper  we  explore  a  decentralized  approach  to  path  generation  for  a  group  of  combat  vehicles  in 
a  battlefield  scenario.  The  mission  is  to  maneuver  the  vehicles  to  cover  a  target  area.  The  vehicles 
are  required  to  maintain  good  overall  area  coverage,  and  avoid  obstacles  and  threats  during  the 
maneuvering.  The  gradient  descent  method  is  used,  where  each  vehicle  makes  its  moving  decision 
by  minimizing  a  potential  function  that  encodes  information  about  its  neighbours,  obstacles,  threats 
and  the  target.  We  conduct  analysis  of  vehicle  behaviors  by  studying  the  vector  field  induced  by  the 
potential  function.  Simulation  has  shown  that  this  approach  leads  to  interesting  emergent  behaviors, 
and  the  behaviors  can  be  varied  by  adjusting  the  weighting  coefficients  of  different  potential  function 
terms. 


1  Introduction 


Autonomous  unmanned  vehicles  (AUVs)  are  receiving  tremendous  interest  due  to  their  potentially  revo¬ 
lutionizing  applications  in  defense,  transportation,  weather  forecast,  and  planetary  exploration  [1],  These 
vehicles  are  often  deployed  in  groups  to  perform  complicated  missions.  Communication  is  often  limited 
in  these  applications  due  to  the  large  number  of  vehicles  involved,  limited  battery  power,  and  constraints 
imposed  by  environmental  conditions  or  mission  requirements.  Hence  a  decentralized  approach  to  co- 
*This  research  was  supported  by  the  Army  Research  Office  under  the  ODDR&E  MURI01  Program  Grant  No.  DAAD19- 
01-1-0465  to  the  Center  for  Networked  Communicating  Control  Systems  (through  Boston  University). 
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ordination  and  control  of  multi-vehicles  is  especially  appealing.  A  decentralized  method  has  another 
advantage  over  a  centralized  one:  it  is  more  robust  to  the  problem  of  single-point  failure. 

Inspired  by  the  emergent  behaviors  demonstrated  by  swarms  of  bacteria,  insects,  and  animals,  control 
methods  which  yield  desired  collective  behaviors  based  on  simple  local  interactions  are  of  great  interest 
[2,  3,  4].  Artificial  potentials  or  digital  pheromones  are  typically  involved  in  such  methods  for  multi¬ 
vehicle  control,  see  e.g.,  [3,  5,  6,  4]  and  the  references  therein.  The  potential  function  method  has  been 
used  in  various  robotic  applications  [7].  The  idea  is  to  derive  a  force  or  other  input  (e.g., velocity)  from 
some  potential  function  which  encodes  relevant  information  about  the  environment  and  the  mission. 

In  this  paper  we  explore  a  decentralized  approach  to  path  generation  for  a  group  of  combat  vehicles 
in  a  battlefield  scenario  using  the  potential  function  method.  The  mission  is  to  maneuver  the  vehicles 
to  cover  a  target  area.  The  vehicles  are  also  expected  to  maintain  good  overall  area  coverage  during  the 
maneuvering,  and  avoid  obstacles  and  threats.  At  every  time  instant  each  vehicle  evaluates  its  potential 
function  profile  and  decides  its  velocity  based  on  the  gradient  descent  method.  The  potential  function 
consists  of  several  terms  reflecting  the  objectives  and  the  constraints.  It  is  constructed  in  such  a  way 
that  only  information  about  neighbouring  vehicles,  local  information  about  dynamic  threats,  and  some 
static  information  (about  stationary  threats,  targets)  are  involved. 

We  analyze  vehicle  formations  at  equilibria  by  taking  the  sum  of  potential  functions  of  all  vehicles 
as  a  Lyapunov  function  candidate.  We  also  study  the  behavior  of  a  vehicle  experiencing  both  attraction 
from  the  target  and  repulsion  from  the  obstacles  by  analyzing  its  vector  field.  Simulation  is  performed 
in  Matlab,  and  it  shows  that  the  decentralized  approach  leads  to  interesting  emergent  behaviors,  and 
the  behaviors  can  be  varied  by  adjusting  the  weighting  coefficients  of  diffrent  potential  function  terms. 

The  remainder  of  the  paper  is  organized  as  follows.  In  Section  2  we  describe  the  problem  setup  and 
construct  the  potential  functions.  We  perform  analysis  of  vehicle  behaviors  in  Section  3.  Simulation 
results  are  reported  in  Section  4.  Section  5  concludes  the  paper. 


2  Potential  Functions 


We  study  the  kinematic  planning  problem  for  N  vehicles  moving  on  a  (two  dimensional)  plane.  Extension 
to  three  dimensional  space  is  straightforward,  although  the  analysis  will  be  more  complicated.  Each 
vehicle  is  treated  as  a  point.  The  coordinates  of  the  i-th  vehicle  Vi  is  denoted  as  pi  =  (xi,  yi),  1  <  i  <  N. 

The  task  for  the  vehicles  is  to  move  toward  and  then  occupy  a  connected  target  area  A  C  M2.  They 
should  avoid  to  crash  into  obstacles  that  are  distributed  in  the  battlefield.  There  are  also  threats,  both 
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stationary  ones  and  moving  ones,  that  endanger  the  vehicles  if  they  are  close.  We  assume  that  each 
vehicle  has  the  knowledge  of  locations  of  stationary  threats,  and  it  can  detect  a  moving  threat  if  the 
threat  is  within  the  distance  Rm. 


The  vehicles  can  talk  to  each  other  and  exchange  information  about  their  positions  and  velocities  if 
they  are  within  the  neighbouring  distance  Rc.  Let  V(t)  be  the  set  of  vehicles  alive  at  t.  For  V)  £  V(t), 
we  define  its  neighbouring  set 

A r(Vi)  =  {V.j  £  V(i)  : ||  Pi  -  Pj  ||<  Rc}. 

If  two  or  more  vehicles  get  too  close,  there  is  a  chance  of  collision.  That  also  makes  it  easier  for  the 
enemy  fire  to  target  the  vehicles.  Another  disadvantage  of  being  too  close  is  that  the  overall  coverage 
area  is  small.  On  the  other  hand,  if  vehicles  are  too  far  apart,  they  lose  contact.  Therefore,  there  is  an 
optimal  distance  ro  <  Rc  between  two  vehicles. 


From  the  above  discussions,  there  are  multiple  objectives/constraints  when  a  vehicle  makes  the 
moving  decision.  To  accomodate  that  we  construct  a  potential  function  J\  for  each  vehicle  Vi  at  time  t, 
where  J}  consists  of  several  terms,  each  term  reflecting  a  goal  or  a  constraint.  To  be  specific, 

Ji  =  XgJ9(Pi(t))  +  XnJn(Pi(t),{Pj(t)}m)  +  X0J°(pi(t))  +  XsJs(Pi(t))  +  XmJm(Pi(t),t),  (1) 


where  J9 ,  Jn,  J°,  Js,  Jm  are  components  of  the  potential  function  relating  to  the  target,  neighbouring 
vehicles,  obstacles,  stationary  threats,  and  moving  threats,  respectively,  and  Xg,  Xn,  X0,  Xs,  and  Xm  >  0 
are  weighting  coefficients.  The  velocity  of  the  vehicle  Vi  is  specified  by 


dJ[ 

dpi ' 

We  now  describe  in  detail  the  components  of  the  potential  function. 


(2) 


•  The  target  potential  J9.  Denote  p(p,A)  =  infa£y4  ||  p  —  a  ||,  the  distance  from  the  point  p  to  the 
target  area  A.  We  then  let 

J9{Pi)  =  f9{p{Pi,A)), 

where  f9{-)  is  a  strictly  increasing  function,  and  f9( 0)  =  0.  This  guarantees  that  in  the  absence  of 
other  objects,  the  vehicle  will  move  toward  the  target.  For  analysis  and  simulation  in  this  paper, 
we  choose  f9(r)  =  r2. 

•  The  neighbouring  potential  Jn.  Since  distance  (as  opposed  to  direction)  is  our  concern  here,  we 
let 

Jn(pi{t),{vmrn)=  E  /"(II K(t)-Pi(t)  II).  (3) 


3 


15 

a 

o  10 


Neighbouring  distance 


Neighbouring  distance 


Figure  1:  An  example  of  the  neighbouring  potential  function. 

where  fn  :  M+  — >  R  is  a  differentiable  function  that  has  the  following  properties:  a)  fn(r )  ap¬ 
proaches  infinity  as  r  — >  0,  and  is  strictly  deceasing  in  [0,ro];  b)  it  is  strictly  increasing  in  [ro,i?c] 
and  ^r(Rc)  =  0.  These  properties  enable  two  vehicles  to  keep  the  optimal  distance  in  the  absence 
of  other  objects,  and  make  the  dynamics  transition  seamless  when  the  neighbouring  set  of  a  vehicle 
is  changing.  An  example  of  such  fn  and  its  derivative  is  shown  in  Figure  1,  where  ro  =  2,  Rc  =  8. 

•  The  obstacle  potential  J°.  An  obstacle  is  a  connected,  closed  set  (could  be  a  single  point)  that  a 
vehicle  cannot  enter.  We  assume  that  there  are  a  finite  number  of  obstacles  {Oj Then 

No 

•Hk)  =  £/>(»,  Oj)),  <4> 

3= 1 

where  p(pi,Oj)  is  the  distance  from  pi  to  the  set  Oj,  and  f°{-)  :  R+  — *  R  is  a  strictly  decreasing 
function  and  f°(r)  — >  oo  as  r  — »  0.  One  example  of  f°  is  f°(r)  =  \. 

•  The  potential  Js  due  to  stationary  threats.  Stationary  threats  can  be  modeled  similarly  as  ob¬ 
stacles,  so  that  vehicles  will  avoid  to  get  close  to  them.  Anisotropic  threats  (dangers  that  are 
direction-dependent)  can  be  taken  care  of  using  appropriate  potential  functions. 

•  The  potential  Jm  due  to  moving  threats.  A  moving  threat  is  a  moving  point  mass.  The  i- th  vehicle 
is  able  to  see  the  moving  threat  Mj  if  ||  Pi  —  qj  ||<  Rd,  and  is  killed  by  Mj  if  ||  pt  —  qj  ||<  Re  <  Rd, 
where  qj  denotes  the  position  of  Mj.  Let  Mi(t)  be  the  set  of  moving  threats  in  the  i-th  vehicles’s 
detection  range,  then  we  let 

Jm(Pi,t)=  Y.  /m(ll  Pi  —  1j  II).  P) 

j:Mj&Mi(t ) 

where  the  function  fm  :  (Re,  oo)  — >•  M  is  differentiable,  strictly  decreasing  on  ( Re,Rd ),  constant  on 
(Rd,  oo),  and  fm(r )  — ►  oo  when  r  — >  Re.  One  can  see  that  with  this  potential  function,  a  vehicle 
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tries  to  keep  at  least  a  distance  Re  from  moving  threats,  and  its  vector  field  remains  continuous 
when  moving  threats  enter  or  leave  its  detection  range.  A  simple  example  for  such  /m(-)  is 


/'»  = 


1  it  r>  ^  ^  Re+Rd 

(r-Re)^  It  Me  <  V  S  2 


16  (r-Rd)2 


8  Re 


(Rd-ReYHRd  +  Re)  ( Rd-Re )S 


if  <  r  <  Rd 


8  Re 


(Rd-Re) 


X  if  r  >  Rd 


3  Qualitative  Analysis  of  Vehicle  Behaviors 


In  this  section,  we  analyze  vehicle  behaviors  under  the  gradient  descent  method.  In  particular,  we 
analyze  how  vehicles  settle  down  after  they  enter  the  target  area,  and  study  the  behavior  of  a  vehicle 
when  it  experiences  both  attraction  from  the  target  and  the  repulsion  from  the  obstacles. 


3.1  Stability  of  equilibrium  configurations 

We  first  consider  multi-vehicles  inside  the  target  area.  We  are  interested  in  knowing  whether  the  vehicles 
will  settle  down  for  an  equilibrium  configuration  under  interactions,  and  if  so,  whether  the  equilibrium 
configuration  is  stable.  Here  we  assume  that  the  only  component  of  the  potential  function  of  each  vehicle 
is  Jn. 


Proposition  3.1  Let  N  be  the  number  of  vehicles. 


1.  For  any  N,  the  configuration  of  vehicles  converges  to  an  equilibrium  under  interactions. 

2.  For  N  =  2,  the  vehicles  maintain  a  distance  of  ro  in  the  equilibrium  configuration  and  the  equilib¬ 
rium  is  globally  stable. 


3.  For  N  =  3,  assuming  that  is  strictly  increasing  in  (0,  ro] ,  there  are  two  possible  equilibrium 
configurations,  equilateral  triangular  with  spacing  ro  (Figure  2  (a)),  and  collinear  with  equal  spacing 
r'  ( Figure  2  (b)),  where  <  r'  <  ro  and 


(6) 


If  is  strictly  increasing  in  [ro,2ro],  r'  is  unique.  The  collinear  configuration  is  unstable,  while 
the  equilateral  configuration  is  stable. 
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Figure  2:  Equilibrium  configurations  for  N  =  3. 


Proof.  1.  Define  a  candidate  Lyapunov  function 


N 


Z=1 


It’s  easy  to  verify  that 


iV 


dJ  „  <9JV  r  .  ,  ll2 

3T  =  -E"  r"'  'll2 


dt  4-f  dp. 
2=1 


(7) 


Hence  J  is  nonincreasing  with  t.  Since  J  is  lower  bounded,  ^  — ►  0.  This  implies  that  pj(t) 
hence  the  vehicles  converge  to  an  equilibrium  configuration. 


0,  Vi,  and 


2.  When  N  =  2,  J(pi(t),p2(t))  =  fn(  ||  pi{t)  -  p2(t)  ||).  Writing  r12(t)  =||  pi{t)  -  p2(t)  ||,  we  have 

dJ  PP ,  , 

Hence  ri2(t)  — >  ro  as  t  — >  oo,  and  it’s  clear  that  this  is  astable  configuration  (as  long  as  ||  Pi(0)—  P2(0)  ||< 

Rc). 

3.  When  N  =  3, 

J(pi(t),p2(t),p3(t))  =  fn(r12(t))  +  fn(r23(t))  +  P(r31(t)), 
where  r^(t)  =||  pi(t)  -  Pj(t)  ||.  Then  ^  = 

„  dp  dp  2  dr  dp  2  dr  r  dfn  2 

-  II  ^r(ri2)ri2  +  -^r(r3i)ri3)  ||  -  ||  -^r(ri2)r2i  +  ^r(r23)r23)  II  -  II  ^r(r3i)r3i  +  ^r(r23)r32)  ||  , 

(8) 

where  f  jj  denotes  the  unit  vector  pointing  from  pj  to  p,.  Eq.  (8)  implies  J  is  strictly  decreasing  unless 
either 


•  ^r(ri2)  =  ^-(r23)  =  ^r(r3i)  =  0,  which  corresponds  to  the  equilateral  triangular  configuration 
shown  in  Figure  2  (a),  or 

•  rjj’s  are  parallel  or  antiparallel  and,  assuming  p2  is  between  p\  and  p3  (we  do  not  lose  generality 
since  three  vehicles  have  symmetric  roles),  ri2  =  r23  =  r'  for  r'  6  (^,ro)  satisfying  (6),  which 
corresponds  to  the  collinear  configuration  in  Figure  2(b). 
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y 

Target 

II'  - p*X 


Obstacle  1 
(-a,-b) 


Obstacle  2 

(a,-b) 
o  Vehicle 


Figure  3:  The  setup  of  two  obstacles  and  one  target 

The  equilateral  triangular  configuration  is  stable.  Indeed,  for  any  perturbation  with  the  corre¬ 
sponding  J  E  (3/n(ro), 2fn(r')  +  fn(2r')),  the  system  will  come  back  to  this  configuration.  The 
collinear  configuration  is  not  unstable,  since  small  perpendicular  perturbation  of  the  middle  vehicle 
leads  to  strict  decrease  of  J  and  the  system  will  converge  to  the  equilateral  triangular  configuration. 


□ 


Remark  3.1  We  note  that  similar  results  appeared  in  [6]  where  the  second  order  dynamics  of  point 
masses  was  considered. 


3.2  Vector  field  analysis 


First  we  consider  the  scenario  as  shown  in  Figure  3.  The  target  is  located  at  the  origin  (0,0).  There 
are  two  (point)  obstacles  located  symmetrically  about  the  y  axis  with  coordinates  (—a,  — b )  and  (a,  — b ), 
respectively,  where  a,  b  >  0.  The  potential  function  in  terms  of  {x,  y )  is  taken  to  be 


As(aE  +  y2)  + 


+ 


(x  +  a)2  +  (y  +  b)2  (x  -  a)2  +  (y  +  b)2  ’ 


and  the  associated  vector  field  is 


x(t) 

y{t) 


2(x+a) _ i _ 2{x-a) _ n\  _ 

[{x+a)^+{y+b)^  [(x-a)^+(y+b)^Y  9 

2fa+b)  | _ 2(y+b) 

[(x+a)2  +  (j/+b)2]2  '  [{x— a)2+(y+6)2]2 


(9) 


We  consider  a  vehicle  on  the  x  axis,  and  study  whether  it  will  move  toward  the  target  under  the 
vector  field  (9)  when  y  <  0  (the  case  y  >  0  is  simpler  and  can  be  studied  similarly).  Due  to  the  symmetry, 
x  =  0,  so  the  real  question  is  whether  y  >  0.  When  x  =  0, 


4  (y  +  b ) 

[(a2  +  {y  +  6)2]2 


2A  gy. 


(10) 
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y 


Target 

Obstacle  1 

Obstacle  2 

(-a,-b)  1 

1  ' 
)' 

(a,-b) 

’  T* 

y  * 

Figure  4:  Vector  field  on  the  tc-axis  for  the  case  of  two  obstacles  and  one  target. 

Let  y  =  y  +  b.  Obviously,  if  y  >  0,  y  >  0.  In  the  following  we  study  the  case  y  <  0,  i.e. ,  y  <  —b. 


Proposition  3.2  There  is  a  unique  solution  y*  £  (--^=,0)  to 

4 y3  -  3 by2  +  a2b  =  0. 


(11) 


Let 


Let  y*  =  y*  —  b.  Then 


_ _ (-i  _ 

9  (a2  +  r2)21 


4  y 


*2 


a2  +  y 


*2 


)• 


•  If  Xg  >  X* ,  y  >  0,  Vr/  <  -b; 

•  If  Xg  =  X* ,  y  >  0  for  y  £  (— oo,  —6)  except  at  y*  where  y  =  0; 


•  If  Xg  <  X* ,  there  exist  y\ ,  y\  such  that  y%  <  y*  <  y\,  and 


V>  0,  ifye  (-00,2/2) 
y  <  0,  if  ye  (y*2,yt) 
y>  0,  if  ye  (y$,-b) 
y  =  0,  if  y  =  y*  or  y*2 


as  illustrated  in  Figure  f.  Furthermore,  y\  (y2,  resp.)  increases  (decreases,  resp.)  with  X  and 
yl  — ►  —6,  y2  — >  —00  as  A  — »•  00. 


Proof.  Let 


Kv) 


4  y 

(a2  +  y2)2 ' 
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Figure  5:  Vector  field  analysis  for  the  case  of  one  obstacle  and  one  target,  (a)  x-component;  (b)  y- 
component;  (c)  total  vector  field. 


Since 

dh  4 (a2  -  3 y2) 

dy  (a2  +  y2)2  ’  1  J 

h{y)  is  strictly  decreasing  on  (—00,—^=),  and  strictly  increasing  on  (— -^=,0).  From  (12),  is  also 
strictly  increasing  on  (— -^=,0).  Graphical  analysis  reveals  that  there  exists  a  unique  A*,  such  that 
the  line  l(y)  =  2A *(y  —  b)  is  tangent  to  the  curve  h(y)  at  a  unique  y*  G  (— -^=,0).  After  algebraic 
manipulations,  one  can  show  that  y*  satisfies  (11)  and  A*  is  defined  by  (3.2).  The  remaining  claims  of 
the  proposition  follow  from  the  graphical  analysis.  □ 


Next  we  investigate  the  motion  of  a  vehicle  in  the  presence  of  one  target  point  (0,0)  and  one  point 
obstacle  (0,  — b ).  Here  no  constraint  on  the  vehicle  position  is  imposed  except  that  we  focus  on  the  region 
y  <  0.  The  vector  field  is 

21  w  “  2 


x  = 


\x2+(y+b )2 


7/  =  2G+b)  _  o\ 

d  [a ;2+(y+6)2]2  ^9  d 

We  will  discuss  x  and  y  separately.  It’s  easy  to  show  that 

x  >  0,  if  (x,  y)  G  {x  <  0  :  x2  +  (y  +  b)2  >  — ?=}  U  {x  >  0  :  x2  +  (y  +  b)2  <  — =} 

V  X9  V  X9 

x  <  0,  if  (x,  y)  G  {x  >  0  :  x2  +  (y  +  b)2  >  — )=}  U  {x  <  0  :  x2  +  (y  +  b)2  <  ^=} 

V  *9 

x  =  0,  if  x  =  0  or  x2  +  (y  +  b)2  =  —j= 
as  shown  in  Figure  5(a).  We  denote  by  C  the  circle  with  radius  — }=  centered  at  (0,  —b). 

For  y,  it’s  straightforward  to  verify 


(13) 


y  >  0  if  x2  +  (y  +  b)2  > 


or  y  >  —b. 
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However,  the  analysis  is  much  more  involved  when  y  <  —b  and  (x,  y)  is  inside  the  circle  C .  The  technique 
used  in  the  proof  of  Proposition  3.2  turns  out  to  be  relevant  here.  Again  let  y  =  y  +  b.  Writing 

2  y 


for  fixed  x.  we  have 


hx{-y)  [x2  +  y2]2  ’ 
dhx  2(x2  -  3 y2) 


dy  (x2  +  y2)3  ' 

From  (14),  we  find  that  hx(y)  is  strictly  decreasing  on  (— oo,  —  ^|x|)  and  strictly  increasing  on  (— ^|x|,  0) 
with  its  derivative  strictly  increasing  on  (— A^|x|,0).  Hence  there  exists  a  unique  A x,  such  that  the  line 
l(y)  =  2A g(y  —  b)  is  tangent  to  the  curve  hx(y)  at  a  unique  yx  €  (— ^|x|,0).  It’s  easy  to  verify  that  yx 
is  a  solution  to 

4  y3  —  3  by2  +  bx 2  =  0,  (15) 


and 


A*  = 


9  16  {{yx)2  —  byx )3 


(16) 


Lemma  3.1  Both  yx  and  Xx  strictly  decrease  as  |x|  increases. 


Proof.  By  the  implicit  function  theorem,  yx  is  a  differentiable  function  of  x2.  Differentiating  (15)  with 
respect  to  z  =  x2  leads  to 

d.vx  -h 

<0, 


dz  12  (yx)2  +  6  byx 
since  yx  <  0.  From  (16),  Xx  also  strictly  decreases  as  \x\  increases.  □ 

From  Lemma  3.1,  for  fixed  Xq  >  0,  there  exists  x\  >  0,  such  that 

A x  =  X g  if  |x|  =  x\ 

X g  >  X g  if  |x|  <  x\  ■ 

Xx  <  X g  if  |x|  >  x\ 


Similarly  as  in  Proposition  3.2,  we  have 


if  A g  >  Xx  (i.e.,  |x|  >  x\),  y  >  0,  Vy  <  - b ; 

if  A g  =  X x  (i.e.,  |x|  =  x\),  y  >  0,  Vy  G  (— oo,  —b)  except  at  yx  =  yXx  —  b  where  y  =  0; 

if  A  g  <  X x  (i.e.,  |x|  <  x\),  there  exist  yf  and  yf  such  that  y^  <yx  —  b  <  yf,  and 

/ 

y  >  0,  if  y  e  (— oo,  yf) 
v  <  o,  if'  y  e  (y^yf) 

y  >  0,  if  y  G  (yf ,  —6) 

V  =  0,  if  y  =  yf  or  yf 


10 


Lemma  3.2  x\  +  ( yXx )2  <  j-.  When 
decreases,  lim|I|_>0y®  =  —b,  and  y®  +  b  > 


x\  <  x\,  yf  (yf, 

1 

\As 


resp.) 


increases  ( decreases,  resp.) 


as  x 


Proof.  As  is  easy  to  verify,  y  >  0  at  (. x ,  y)  if  x2  +  (y  +  b )2  >  j-.  Since  y  =  0  at  (x\,  yXx ),  the  first  claim 
of  the  lemma  holds. 


Since  for  any  y  <  0,  hXl(y)  >  hX2(y )  if  |aq|  >  \x2\,  graphical  analysis  reveals  that  yf  (yf,  resp.) 
increases  (decreases,  resp.)  as  |x|  decreases,  and  limui^o  2/f  =  —  b.  Again,  since  y  >  0  at  (. x,y )  if 
x2  +  (y  +  b)2  >  j-,  the  point  (0,  y$)  is  inside  the  circle  C,  i.e.,  y®  +  b  > - 4=.  □ 

9  V  ^9 


Figure  5  (b)  and  (c)  sketch  the  y-conrponent  of  the  vector  field  and  the  total  vector  field,  respectively. 
We  can  see  that  the  only  point  where  x  =  y  =  0  is  (0,  y®).  But  this  is  an  unstable  equilibrium  as  one  can 
tell  from  Figure  5(c).  We  can  also  verify  that  the  linearized  system  at  (0,  y®)  has  a  positive  eigenvalue. 


4  Simulation  Results 


In  this  section  we  present  some  results  of  simulation  performed  in  Matlab.  As  shown  in  Figure  6(a), 
we  have  ten  vehicles  (the  dots),  two  obstacles  (the  solid  circles),  one  target  (the  dashed  circle),  and  one 
moving  threat  (represented  by  the  cross).  The  moving  threat  is  assumed  to  move  around  the  target 
with  a  constant  angular  velocity  u).  We  don’t  assume  any  stationary  threats  in  the  simulation.  The 
velocity  magnitude  is  bounded  by  1.  The  Matlab  function  “fconmin”  is  called  to  solve  the  constrained 
optimization  problem  for  each  vehicle  at  every  time  instant. 

Figure  6  and  Figure  7  show  snapshots  of  the  trajectories  of  the  vehicles.  The  weighting  constants 
used  in  Figure  6  are  \g  =  50,  Xn  =  5,  A0  =  105,  \m  =  104.  The  angular  velocity  of  the  moving  threat  is 
to  =  10.  We  observe  that  due  to  the  relatively  high  weight  on  obstacles,  the  vehicles  take  long  paths  to 
avoid  vehicles;  however,  if  we  change  A0  to  104,  they  take  the  shorter  path  passing  the  “potential  valley” 
between  two  obstacles  (Figure  7). 

Figure  8  shows  the  effects  of  Am.  When  Am  =  104,  vehicles  can  successfully  enter  the  target  area 
(Figure  8  (a)),  while  they  fail  to  do  so  when  Am  =  105  (Figure  8  (b)). 

From  the  simulation  results,  we  see  that  the  decentralized  approach  based  on  potential  functions 
lead  to  some  emergent  behaviors  of  vehicles.  In  addition,  we  can  modify  the  behaviors  by  appropriately 
changing  the  weighting  constants. 
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Figure  6:  Snapshots  for  the  case  A0 


105:  (a)  Initial  positions;  (b)  Vehicles  detouring  to  avoid  obstacles. 


Figure  7:  Snapshots  for  the  case  A0  =  104:  (a)  Initial  positions;  (b)  Vehicles  passing  the  “potential 
valley”  between  obstacles. 
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x-axis 


x-axis 


(a) 


(b) 


Figure  8:  Effect  of  the  parameter  Xm:  (a)  Xm  =  104  (vehicles  enter  the  target  area);  (b)  Xm  =  105 
(vehicles  fail  to  enter  the  target  area). 


5  Conclusions 


In  this  paper  we  have  developed  a  decentralized  approach  to  coordination  and  control  of  multi-vehicles 
using  potential  functions.  We  designed  potential  function  terms  for  different  objectives  or  constraints.  A 
battlefield  mission  scenario  was  considered,  where  the  vehicles  were  required  to  occupy  a  target  area  (or 
point),  avoid  obstacles,  evade  threats,  and  maintain  reasonable  inter- vehicle  distances.  Preliminary  anal¬ 
ysis  of  vehicle  behaviors  was  presented.  Simulation  was  conducted  and  interesting  emergent  behaviors 
were  observed. 

The  most  important  advantage  of  this  approach  is  simplicity  since  only  local  and  static  information 
is  needed  in  the  path  generation.  It  is  also  flexible  and  robust,  which  is  of  vital  importance  in  complex, 
dynamic  environments  such  as  the  battlefields.  The  disadvantage  is  that  the  possibility  of  being  trapped 
in  local  minima  exists,  which  has  been  a  long  time  concern  in  the  studies  of  the  potential  function  method 

[8] .  Practically  interactions  between  vehicles  and  dynamic  changes  in  the  environment  may  prevent  a 
vehicle  from  being  trapped.  Artificially  introduced  perturbation  will  also  help  to  resolve  this  problem 

[9] , 


Ongoing  work  includes  analysis  of  vehicle  behaviors  in  the  presence  of  vehicle  interactions  as  well  as 
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attractions/repulsions  from  other  objects. 
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